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Abstract 

Close-packed  array  (Bi,Sb)2(Te,Se)3  compound  semiconductor  thermoelectric  power  modules  have  been  used  extensively  in 
terrestrial  thermoelectric  generators  for  converting  thermal  energy  directly  into  electricity.  In  this  paper  a  diffusion  controlled  degrada¬ 
tion  model  based  on  a  new  modified  Whipple-type  grain  boundary  diffusion  model  has  been  developed.  This  model  considers  the 
effect  of  temperature  gradient  on  a  high  diffusivity  path  and  has  been  developed  for  studying  the  degradation  mechanism  of 
(Bi,Sb)2(Te,Se)3  compound  semiconductor  power  thermoelectric  modules  caused  by  the  formation  of  dark  bands. 
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1.  Introduction 

The  Kadel  thermoelectric  energy  conversion  device 
generally  operates  at  hot  junction  temperature  of 
260-316  °C,  and  maintains  temperature  gradients  of 
170°C  between  the  hot  and  cold  junctions.  Owing  to 
the  coefficient  of  thermal  expansion  mismatch  of 
different  parts  of  the  module,  thermal  stress  can  occur 
at  the  interfaces  in  the  thermoelectric  module  system 
during  manufacturing  and  operation.  Surface  morpho¬ 
logy  and  energy  dispersion  X-ray  spectroscopy  studies 
have  shown  that  edge  microcracks,  adjacent  to  the 
insulator  in  the  semiconductor,  propagate  along  the 
interphase  boundary  between  the  oxidized  and  non- 
oxidized  areas.  This  propagation  is  a  result  of  stress 
induced  impurity  diffusion  from  the  microcrack  tip  to 
the  interphase  boundary  ahead  of  the  tip.  This  impurity 
stress-assisted  diffusion  cracking  involves  structural 
changes  that  weaken  and  break  the  atomic  bonds  at  the 
interphase  boundary  zone  ahead  of  the  crack  tip. 
Recently,  a  new  dynamic  stress  assisted  diffusion 
cracking  model  has  been  developed  for  describing  the 
influence  of  thermal  stress  on  the  diffusion  rate  of  the 
impurity  (Fe)  ahead  of  the  crack  tip  and  microcrack 
growth  in  the  brittle  (Bi,  Sb)2(Te,Se)3  semiconductor 
thermoelectric  materials  [1].  Moreover,  electrical 


resistance  and  surface  morphology  studies  further 
indicate  that  the  high  electrical  resistance  is  strongly 
related  to  the  presence  of  an  oxidized  area  of  semi¬ 
conductor  near  the  hot  junction,  which  was  operated 
above  260  °C,  as  well  as  the  dark  band  along  the  inter¬ 
phase  boundary  between  the  oxidized  and  non- 
oxidized  areas  of  the  P  type  semiconductor 
thermoelements.  The  formation  of  the  dark  band  is 
caused  by  surface  diffusion  of  impurity  atoms  from  the 
insulator  (separator)  into  the  microcrack  which  pro¬ 
pagate  along  the  interphase  boundary  inside  the  semi¬ 
conductor.  A  schematic  diagram  and  scanning  electron 
micrograph  showing  the  dark  band  are  illustrated  in 
Fig.  1.  The  major  goal  of  this  study  is  to  investigate  the 
mechanism  of  diffusion  controlled  dark  band  forma¬ 
tion  on  a  common  phenomenological  basis.  Mostly, 
surface  diffusion  of  impurities  along  the  microcrack 
(such  as  Si  with  low  diffusivity)  can  be  logically 
modeled  by  the  same  approach  used  for  grain 
boundary  diffusion.  Accordingly,  the  formation  of  the 
dark  band  due  to  the  surface  diffusion  along  the 
microcrack  high  diffusivity  path  can  be  successfully 
modeled  by  a  new  modified  Whipple-type  grain 
boundary  diffusion  solution  which  takes  into  account 
the  effect  of  temperature  gradient  driving  force  on 
grain  boundary  diffusion.  A  geometric  model  for  grain 
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Fig.  1.  Kadel  thermal  electric  module:  (a)  schematic  diagram  of 
the  power  thermoelectric  module  and  generator  and  (b)  scanning 
electron  micrograph  showing  the  dark  band  in  the  semiconduc¬ 
tor  thermoelement. 


boundary  diffusion  in  the  thermoelectric  module  with  a 
microcrack,  which  is  treated  as  a  grain  boundary,  is 
shown  in  Fig.  2. 
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2.  Physical  model  and  equations 


(b) 


Based  on  the  linear  diffusion  theory,  the  flux  with  a 
solid  state  thermomigration  driving  force  can  be 
expressed  [2]  as 


Fig.  2.  Schematic  illustration  of  (a)  location  of  microcrack  in  the 
Kadel  thermal  electric  module  and  (b)  geometric  model  for  grain 
boundary  diffusion  in  the  thermoelectric  module  with  the 
microcrack  as  a  high  diffusivity  path. 


J= -D\VC- 


(1) 


where  £>,  C,  F  and  k  represent  the  diffusivity,  concen¬ 
tration,  thermomigration  driving  force  and  Boltz¬ 
mann's  constant  respectively. 

Finite  element  calculation  has  shown  that  the 
temperature  profile  for  (Bi,  Sb)2(Te,Se)3  semiconductor 
thermoelectric  power  modules  can  be  treated  reason¬ 
ably  as  linear  in  the  dark  band  [3].  Therefore,  the 
driving  force  for  thermal  diffusion  can  be  written  [4] 


Q^dT_  _Q*AT 
T  dx~  T  Ax 


(2) 


where  Q*  is  the  heat  of  transport  and  dT/dx  is  the 
temperature  gradient.  Then  the  diffusion  with  the 
driving  forces  becomes 


dC  ,  DF 

—  =  DV‘  C - VC 

dt  KT 


(3) 
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Our  approach  to  the  diffusion  problem  is  similar  to 
that  used  by  Whipple  [5].  In  this  paper,  a  two-dimen¬ 
sional  grain  boundary  diffusion  model  influenced  by  a 
thermomigration  driving  force  is  taken  into  account. 
This  driving  force  is  perpendicular  to  the  microcrack. 
Let  the  concentration  and  diffusivity  within  the  high 
diffusivity  path  by  given  by  C2(x,y,t)  and  D2  and  the 
drain  into  the  adjacent  lattice  by  Cx[x,y,t)  and  D,.  The 
width  of  the  high  diffusivity  path  is  <5.  D,  and  D2  are 
logically  assumed  to  be  constant  and  the  temperature 
gradient  independent  within  the  very  narrow  dark 
band.  The  diffusion  equations  with  the  perpendicular 
thermal  gradient  driving  force  can  be  expressed  as 


(tfcy  a^c,\ 

D  yF  1  0C, 

_dc, 

W  dy2l 

KT  dx 

dt 

l&Ci  tfCi) 

D2F2  dC2 

_dC2 

2\  dx2  dy2  1 

KT  dx 

dt 

and  can  be  solved  by  means  of  Fourier-Laplace  trans¬ 
form. 

Consequently  an  exact  solution,  taking  into  account 
the  effect  of  thermal  gradient  driving  force,  can  be 
written 

C{x,y,t)  =  C[{x,y,t)  +  Ch(x,y,t) 


C[x,y,t)=  C„  erfc  |  j 
z  2 -in  i 


da 
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with  dimensionless  parameters 


(9) 


where  F,  and  F2  are  the  effective  driving  forces  in  the 
lattice  and  at  high  diffusivity  path,  respectively. 

Assume  that  the  concentration  and  flux  are 
continuous  at  the  interface  between  the  boundary  and 
the  lattice  (x  =  ±  d/2),  and  that  the  surface  concentra¬ 
tion  maintains  constant  value  Cn. 


,  x  ~{6/2) 

§  (ZV)1/2 
V={D^fr2 


OX  ox 


(5) 
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Rearranging  terms  and  ignoring  terms  of  the  order  of 
62,  the  new  governing  equations  become 


0C=  /a2g'  a2g"\  D\F\  dc 
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where  C  represents  the  total  lattice  concentration 
which  is  due  to  ordinary  diffusion  in  the  lattice  and 
diffusion  out  from  the  high  diffusivity  path  to  the 
lattice.  Eq.  (7)  is  subjected  to  the  initial  and  boundary 
conditions  as  given  by 

C(x,0,0)=C()  C(*,y,0)  =  0 

\o) 

C(x,0,t)=  C0  C(°°,oo,f)  =  0 


A 
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(10) 
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j8  =  (A-l  )a 

where  C,  is  the  lattice  concentration  owing  to  planar 
lattice  diffusion  in  the  absence  of  the  high  diffusivity 
path  contribution,  Ch  is  the  lattice  concentration  owing 
to  out-diffusion  from  the  high  diffusivity  path,  erfc  is 
the  complementary  error  function  and  a  is  the  integra¬ 
tion  variable. 


3.  Discussion 

The  effect  of  the  driving  force  caused  by  the 
temperature  gradient  on  grain  boundary  diffusion 
becomes  evident  by  comparing  Eq.  (9)  with  Whipple’s 
solution.  In  Cx,  the  rj  dependence  in  the  comple¬ 
mentary  error  function  remains  unchanged.  Hence  the 
driving  force  caused  by  temperature  gradient  perpen¬ 
dicular  to  the  high  diffusivity  path  has  no  effect  on  the 
contribution  to  Cj  from  direct  volume  or  lattice  diffu¬ 
sion.  Lattice  diffusion  owing  to  out-diffusion  from  the 
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high  diffusivity  path  Ch  is  explicitly  modified  along  the 
x-direction  parallel  to  the  driving  force  by  the  term 
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in  a  complementary  error  function.  While  atomic  flow 
across  the  high  diffusivity  path  is  more  sensitive  to 
temperature  gradient  than  in  the  lattice,  F2  is  larger 
than  F, .  In  addition,  the  diffusivity  of  impurity  within 
the  high  diffusivity  path  D2  is  much  larger  than  inside 
the  lattice  or  volume  Dx,  so  that  D2FX  is  larger  than 
D2FX  .  Substituting  Eq.  (2)  into  (11),  the  shifting  term  % 
in  the  complementary  error  function  of  Ch  is  given  by 


X 


(7-1 

A  -  1 


D\Qi*  ~  DiQ*  A7-1 
KT2  Ax 


I  (Dru2t112) 


(12) 


where  Q*  is  the  heat  of  transport  in  the  lattice  and  Q* 
is  the  heat  of  transport  in  the  high  diffusivity  path.  Note 
that  Q2*  >  Q*>  0  and  D2  Q*  >  Dx  Q*,  which  indicate 
that  the  shifting  term  %  in  the  argument  of  the  comple¬ 
mentary  error  function  (erfc)  is  negative.  Consequently, 
the  decrease  in  the  argument  of  the  erfc  term  owing  to 
the  negative  shifting  term  x  enhances  impurity  lattice 
penetration  from  the  high  diffusivity  path  into  the  non- 
oxidized  area  below  the  microcracks  treated  as  the 
high  diffusivity  path.  Alternatively,  an  increase  in  the 
argument  of  the  erfc  term  owing  to  the  positive  shifting 
term  %  diminishes  the  impurity  lattice  penetration  from 
the  high  diffusivity  path  into  the  oxidized  area  above 
the  microcrack.  However,  the  overall  lattice  composi¬ 
tion  profile  is  adjusted  by  the  shifting  term  along  the 
direction  of  thermomigration.  It  is  worth  noting  that 
the  formation  of  the  dark  band  caused  by  redistri¬ 
bution  of  impurities  along  the  microcrack  (high  diffu¬ 
sivity  path)  and  in  the  matrix  (lattice)  of  the 
semiconductor  is  influenced  not  only  by  the  tempera¬ 
ture  gradient  (thermomigration),  but  also  by  the 
electric  current  (electromigration)  in  the  direction 
normal  to  the  high  diffusivity  path.  Electromigration 


may  reinforce  or  counteract  formation  of  the  dark 
band,  which  depends  on  the  direction  of  electric 
current  and  the  type  of  semiconductor  (N-  or  P-type) 
involved. 


4.  Conclusions 

This  diffusion  degradation  model  is  in  good  agree¬ 
ment  with  the  dark  band  observed  in  the  brittle 
(Bi,Sb)2(Te,Se)3  thermoelements.  There  is  a  strong 
correlation  between  the  temperature  gradient  driving 
force  and  the  formation  of  a  dark  band  that  degener¬ 
ates  the  thermoelectric  power  and  increases  the  elec¬ 
trical  resistance  of  the  thermoelectric  module.  A 
thermomigration  driving  force  enhances  penetration  of 
impurities  into  semiconductor  thermoelements  and 
accelerates  the  rate  of  performance  degradation.  This 
model  considers  the  effect  of  the  thermomigration 
driving  force  perpendicular  to  high  diffusivity  path  and 
may  also  be  applicable  to  ceramic  materials  with 
temperature  differences  between  the  surface  and 
interior  regions  owing  to  poor  thermal  conductivity. 
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